ON SOME ERRORS RELATED TO THE GRADUATION OF 
MEASURING INSTRUMENTS 



Nicolas Bouleau* 
October 2006 



Abstract. The error on a real quantity Y due to the graduation of the measuring 
instrument may be approximately represented, when the graduation is regular and 
fines down, by a Dirichlet form on K. whose square field operator does not depend 
on the probability law of Y as soon as this law possesses a continuous density. This 
feature is related to the "arbitrary functions principle" (Poincare, Hopf). We give 
extensions of this property to M. d and to the Wiener space for some approximations 
of the Brownian motion and apply these results to the discretization of stochastic 
differential equations. 
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I. Introduction. 

1.1. The approximation of a random variable Y by an other one Y n yields most 
often a Dirichlet form. The framework is general, cf. Bouleau (2006) : Let Y and 
Y n be defined on (fl, A, P) with values in the measurable space (E, J 7 ), denoting Py 
the law of Y, if there exists an algebra T> of bounded functions dense in L 2 (Py) 
and a sequence of positive real numbers a n such that, for all (p £ T> there exists 
AM e L 2 (Py) such that 



Wip G T> a n E[(<p(Y n ) - V {Y)){^{Y n ) - ^{Y))\ - -2 < %>], V 



-) 



then £[ip, ip\ = — < A[tp], ip > is a Dirichlet form. 

Often, when this Dirichlet form exists and does not vanish, the conditional law of 
Y n given Y = y is not reduced to a Dirac mass, and the variance of this conditional 
law yields the square field operator T. On the other hand when the approximation 
is deterministic, i.e. when Y n is a function of Y say Y n = r] n (Y), then most often the 
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symmetric bias operator A and the Dirichlet form vanish, Bouleau (2006, examples 
2.1 to 2.9 and remark 5). 

Nevertheless, there are cases where the conditional law of Y n given Y is a Dirac 
mass, i.e. Y n is a deterministic function of Y, and where the approximation of Y by 
Y n yields even so a non zero Dirichlet form on L 2 (Py). 

This phenomenon is interesting, insofar as randomness (here the Dirichlet form) 
is generated by a deterministic device. In its simplest form, the phenomenon appears 
precisely when a quantity is measured by a graduated instrument when looking for 
the asymptotic limits as the graduation fines down. 

We will first expose (section 1.2) the simplest case of a real quantity measured 
with an equidistant graduation. The mathematical argument is here the same as 
for the arbitrary functions method, about which we give a short historical survey 
(section 1.3). 

The second part is devoted to theoretical tools that we shall use later on. We 
first recall the bias operators and the properties of the Dirichlet form associated 
with an approximation. Next we prove a version of a Girsanov theorem for Dirichlet 
forms which has its own interest, i.e. an answer to the question of an absolutely 
continuous change of measure for Dirichlet forms. At last we recall some simple 
properties of Rajchman measures. 

Thanks to these tools, in the third part we take up again the classical case in 
order to extend it to M. d and to more general graduation. 

The fourth part is concerned by Rajchman martingales, i.e. continuous local 
martingales whose brackets are a.s. Rajchman measures. These martingales possess 
remarkable properties of weak convergence for some approximations. Here the fol- 
lowed method is essentially an extension of the seminal idea of Rootzen (1980). If 
we restrict the framework to the Wiener space, the obtained limit quadratic forms 
may be shown to be closable hence Dirichlet, this is done in the fifth part. This 
yields asymptotic speed of weak convergence for the discretization of stochastic dif- 
ferential equations in the case encountered for mechanical systems which was not 
explicited in the recent works of Kurtz and Protter (1991), Jacod and Protter (1998). 

1.2. The basic example. 

Let Y be a real random variable. It is approximated by Y n to the nearest grad- 
uation, i.e. 

n 2n 

([x] denotes the entire part of x, and {x} — x — [x] the fractional part). 

We put Y n = Y + £ n (F) where the function £ n (x) = ^ — ^ — x is periodic with 
period ^ and may be written £ n (x) = ^9(nx) with 6{x) — \ — {x}. 
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fig 1 and 2: Approximation of the quantity Y 
to the nearest graduation. 



Let Py the law of Y, we approximate Y by Y n on the algebra V = C 1 n Lip with 
the sequence a n = n 2 , cf. Bouleau (2006). 

Let us recall that the Rajchman class is the set of bounded measures on K. whose 
Fourier transform vanishes at infinity. These measures are continuous (do not charge 
points) and are a band in the space of bounded measures on H. (see section 11.3 below). 



Theorem 1. //Py is a Rajchman measure, 

(n(Y n -Y),Y) =4 (V,Y) 
where V is uniform on (— |, \) independent ofY, and for (p E C 1 fl Lip 

1 
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Ey[^]. 



(2) 



Here denotes the weak convergence, i.e. the convergence of probability measures 
on continuous bounded functions, Ey is the expectation under Py. 
Proof. It is equivalent to study the weak convergence of (| + n(Y n — Y), Y) = 
(| + 9(nY), Y). Since | + 6 takes its values in the unit interval, it is enough to study 
the convergence on the characters of the group T 1 x R, i.e. 



E[e 



2mk{\+8{nY)) JuY] 



E[e 



- 2iirknY ^iuY~i 



^y(w - 2-Kkn) 



where \l/y is the characteristic function of Y . This tends to ^/(u)l^o} since Py is 
Rajchman, proving the first assertion. 

If if E C 1 n Lip, the relation <p(Y n ) - <p(Y) = (Y n - Y) ^ <p'(Y + a(Y n -Y))da 
gives 

n 2 E[( V (Y n ) - V {Y)) 2 } = E[9 2 (nY) V ' 2 (Y)] + o(l) 



and E[6 2 (nY)(p' 2 (Y)] -> f\ 9 2 (t)dtE[ip' 2 (Y)) what proves the second one. 



o 
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If Py is absolutely continuous and satisfies the Hamza condition, Fukushima et 
al. (1991, theorem 3.1.6 p. 105), e.g. as soon as Py has a continuous density, the form 
£[<f\ = ^Ey[(p' 2 ] is Dirichlet and admits the square field operator T[ip] = j^p' 2 - The 
graduation yields therefore an error structure (R, B(M), Py, D, T), Bouleau (2003, 
Chapter III) whose operator T does not depend on Py provided that Y has a regu- 
lar density. This translates in terms of errors the arbitrary functions principle. 

Remark. In his study of the exception to the law of errors of Gauss, Henri Poincare 
(1912, p. 219) considers the case where the nearest graduation is chosen except in a 
central zone across the middle between two marks of the graduation, where either 
the left one or the right one is chosen randomly with equal probability. It is easily 
seen that this case may be handled in the same way and gives similar results. The 
Dirichlet form is increased with respect to the above case. 

1.3. Historical comment. 

In his intuitive version, the idea underlying the arbitrary functions method is 
ancient. The historian J. von Plato (1983) dates it back to a book of J. von Kries 
(1886). We find indeed in this philosophical treatise the idea that if a roulette 
had equal and infinitely small black and white cases, then there would be an equal 
probability to fall on a case or on the neighbour one, hence by addition an equal 
probability to fall either on black or on white. But no precise proof was given. The 
idea remains at the common sense level. 

A mathematical argument for the fairness of the roulette and for the equi- 
distribution of other mechanical systems (little planets on the Zodiac) was proposed 
by H. Poincare in his course on probability published in 1912 (1912, Chap. VIII 
§92 and especially §93). In present language, Poincare shows the weak convergence 
of tX + Fmod 2ir when t \ oo to the uniform law on (0, 2ir) when the pair (X, Y) 
has a density. He uses the characteristic functions. His proof supposes the density 
be C 1 with bounded derivative in order to perform an integration by parts, but 
the proof would extend to the general absolutely case if we were using instead the 
Riemann-Lebesgue lemma. 

The question is then developed without major changes by several authors, E. 
Borel (1924) (case of continuous density), M. Frechet (1921) (case of Riemann- 
integrable density), B. Hostinski (1926) (1931) (bidimensional case) and is tackled 
anew by E. Hopf (1934), (1936) and (1937) with the more general point of view of 
asymptotic behaviour of dissipative dynamical systems. Hopf has shown that these 
phenomena may be mathematically understood in the framework of ergodic theory 
and are related to mixing. Today the connection is close to Rajchman (or mixing) 
measures cf. Lyons (1995), Katok and Thouvenot (2005), interesting objects related 
to deep properties of descriptive set theory. 
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II. Theoretical tools 



1 1.1. Approximation, Dirichlet forms and bias operators. 

Our study uses the theoretical framework concerning the bias operators and the 
Dirichlet form generated by an approximation proposed in Bouleau (2006). We 
recall here the definitions and main results for the convenience of the reader. Here, 
considered Dirichlet forms are always symmetric. 

Let Y be a random variable defined on (ft, A, P) with values in a measurable 
space (E, J-) and let Y n be approximations also defined on (Q, A, P) with values 
in (E, J 7 ). We consider an algebra T> of bounded functions from E into R or C 
containing the constants and dense in L 2 (£, JF, Py) and a sequence a n of positive 
numbers. With V and (a n ) we consider the four following assumptions defining the 
four bias operators 

f VifEV, there exists A[<p] G L 2 (E, J 7 , P y ) s.t. V^eP 
\ lim n _a n E[( V (y n ) - <p(Y)) X (Y)] = E y L%] X ]. 

( VcpeV, there exists A[tp\ G L 2 (E, .F,Py) s.t. V^eP 
\ lim n ^ 00 a n E[( V 9(r) - <p(Y n )) X (Y n )] = E Y [A[ip} X }. 

Mip E V, there exists A[(p] E L 2 {E,T, Py) s.t. Vx G V 
lim n _ >0O a n E[(^(y n ) - <p{Y)){ X {Y n ) - X (Y))} = -2E Y [A[<p] X ]. 

V(p G V, there exists A[(p] G L 2 (E, J 7 , P y ) s.t. \/ X G V 
lim^oo a„E[(v9(F n ) - ^y))(x(y„) + xM)] = 2E y ^Mx]- 

We first note that as soon as two of hypotheses (HI) (H2) (H3) (H4) are fulfilled 

(with the same algebra V and the same sequence a n ), the other two follow thanks 

to the relations _ _ 

~ A+A A-A 

2 T 2 

When defined, the operator A which considers the asymptotic error from the point 
of view of the limit model, will be called the theoretical bias operator. 

The operator A which considers the asymptotic error from the point of view of 
the approximating model will be called the practical bias operator. 

Because of the property 

< A[lf], X > L 2 {rY) = < <P,A\X] >L2(Py) 

the operator A will be called the symmetric bias operator. 

The operator )\ which is often (see theorem 3 below) a first order operator will 
be called the singular bias operator. 

Theorem 2. Under the hypothesis (H3), 
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a) the limit 

— f)< 

£[<p, X ] = lim -f E[(<p(Y n ) - <p(Y))(x(Y n ) - X (Y)] (3) 

defines a closable positive bilinear form whose smallest closed extension is denoted 
(5,D). 

b) (£,©) is a Dirichlet form 

c) (£, D) admits a square field operator T satisfying Vy?, x^V 

V\p\=A\^\-%pA^p (4) 
E Y [I>] X ] = \ima n EMY n ) - <p(Y)) 2 ( X (Y n ) + x(Y))/2] (5) 

n 

d) (£, D) is local if and only ifW^p G V 

lima n E[( V (Y n )- V (Y)) 4 ] =0 (6) 

n 

this condition is equivalent to 3A > 2 lim n a n E[|<^(Y^) — y?(F)| A ] = 0. 

e) If the form (£, D) is local, then the principle of asymptotic error calculus is 
valid on V = {F(f f p ) : fceV, F e C 1 (1R P , R)} i.e. 

lim n a n E[(F(/ 1 (y n ), • • • , " F(fi(Y), • • • , / P Q0) 2 ] 

= EyEU Wi, • • • - /p)^f(/i, • • • - /p)r[/i, /;]]■ 

An operator £> from V into L 2 (Py) will be said to be a /irst order operator if it 
satisfies 

B[i PX \ = B[<p\x + vB\ K \ y v ,xev 

Theorem 3. Under (HI) to (H4). If there is a real number p > 1 s.i 
limaJE[(^(y n ) - p(y))XY n ) - if>(Y)\>] = V</?, ip ET> 

n 

then A is first order. 

In particular, if the Dirichlet form is local, by the d) of theorem 2, the operator 
A is first order. 

Example 1. In order to deepen the discussion begun above in the introduction 
about the conditional law of Y n given Y = y, we give below a simple example where 
both the conditional law of Y n given Y = y and the conditional law of Y given 
Y n = y are Dirac measures and where nevertheless the approximation gives rise to 
a non-zero Dirichlet form. 

Let us consider the unit interval and the dyadic representation of real numbers. 
If Y is uniformly distributed we may write Y = YlkLo 2^+1 wnere the are inde- 
pendent identically distributed with law |<5 + \5\. 
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Let us approximate Y by Y n = J2k=o 2^r + \ J2T=n wrr- We see tliat Y and Y n 
are deterministically linked : 

2 2 n 2 2 n 

Now, it is easily seen that on the algebra V = £{e 2tnkx , 6 Z} we have 

3A n E[( V (Y n ) - <p{Y))ty{Y n ) - ${Y))\ -> E[<^], 

what gives in the real domain the Dirichlet form £[(p] = |E[</?' 2 ]. o 

11.2. Girsanov theorem for Dirichlet forms. 

An error structure is a probability space (Q, A, P) equipped with a local Dirichlet 
form with domain D dense in L 2 (fl,A, P) admitting a square field operator T, see 
Bouleau (2003). We denote PA the domain of the associated generator. 

Theorem 4. Let (Q, A, P, O, T) &e an error structure. Let be f G D n L°° snc/i i/iat 
/ > 0, E/ = 1, W^e P x = /.P. 

aj T/ie bilinear form £ 1 defined on VA n L°° fry 



5i[u,u] = — E 
closable in L 2 (Pi) and satisfies for u, v G "DA fl L 



(7) 



1 

£i[w,n] = - < >= - < w, A\v >= -E[fT[u,v}] (8) 

where Ai[u] = A[u] + jjT[u, /]. 

6^) Lei (Di, £1) be the smallest closed extension of (VA fl L°°,£i). Then D C D 1; 
£1 is local and admits a square field operator Ti, and 

T 1 = T on D 

in addition VA C VA 1 and Ai[u] = A[u] + ^jT[u, f] for all u G VA. 

Proof. 1) First, using that the resolvent operators are bounded operators sending 
L°° into VA fl we see that VA fl L°° is dense in D (equipped with the usual 
norm (||.|| 2 2 + £[.]) 1//2 ), hence also dense in L 2 (Pi). 

2) Using that D fl L°° is an algebra, for u, v G VA fl L°° we have 

£i[u,v] = -E[fvA[u] + \vY\uJ}} = h:[T[fv,u}-vT[u,f]] = ±E[fV[u,v]]- 
So, defining A 1 as in the statement, we have Vit, v G ©A fl L°° 

^[ttju] = -Ei[uAiu] = -Ei[itAiu]. 
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The operator A\ is therefore symmetric on VA fl L°° under Pi. Hence the form E\ 
defined on VA fl L°° is closable, Fukushima et al. (1994, condition 1.1.3 p 4). 

3) Let (H>i,£i) be the smallest closed extension of (VA fl L°° ,Si). Let be u G D 
and u n G VA fl with u n — > w in D. Using £i[u n — u m ] < ||/||oo^[ M n — M m] and 
the closedness of £i we get u n — > u in D 1; hence DcDj. Now by usual inequalities 
we see that T[u n ] is a Cauchy sequence in L 1 (Pi) and that the limit Tif-u] does not 
depend on the particular sequence (u n ) satisfying the above condition. Then fol- 
lowing Bouleau (2003, Chap. Ill §2.5 p. 38) the functional calculus extends to Di, 
the axioms of error structures are fulfilled for (fl, A, Pi, Di, Ti) and this gives with 
usual arguments the b) of the statement. o 

11.3. Rajchman measures. 

Definition 1. A measure p on the torus T 1 is said to be Rajchman if 



The set of Rajchman measures 1Z is a band : if p G 1Z and if v <C \p\ then v G 7Z, 
cf. Rajchman (1928) (1929), Lyons (1995). 

Lemma. Let X be a real random variable and let ^x(u) = Ke tuX be its character- 
istic function. Then 




when \n\ j oo. 





Proof, a) If lini| uHoo V x (u) = then V x (2nn) = (F {x} }(n) -> 0. 
b) Let p be a probability measure on T 1 s.t. p G 1Z. From 




oo 



we have 




■ p e 2i7T[u]x p(dx). 



Since x p p(dx) G 7£ 




with \a p \ < 1 and lim| n |_ >0O a p (n) = 0, so 
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Now if P{x} e 72, since l {:re [p,p + i[}.P{x} < P{x} we have 

lim E[e 2M ] = lim £ E[e 2M l {Xe[p , p+1[} ] 

which goes to zero by dominated convergence. o 

A probability measure on R satisfying the conditions of the lemma will be called 
Rajchman. 

Examples. Thanks to the Riemann-Lebesgue lemma, absolutely continuous mea- 
sures are in 72. It follows from the lemma that if a measure v satisfies v ★ • • • * v e 72 
then v e 72. 

Let /3 e]0, |[, let K x , be the closed set obtained by taking away from the unit 
interval the open interval of length 1 — 2(3 centered on 1/2, the closed set ob- 
tained by iterating these suppressions homothetically on each segment of K^-i and 
let /ifc be the probability measure j^l K k -dx The continuous measure \i weak limit 
of the /ik carried by the perfect set C\kK k is in 72 iff 1/(5 is not a Pisot number (a 
Pisot number is a root of a polynomial with entire coefficients and with coefficient 
of highest degree term equal to 1, irreducible over Q, such that the other roots have 
a modulus < 1) cf. Kahane and Salem (1953). 

Proposition 1. Let X, Y, Z be random variables with values in R 7 R, and R m resp. 
Then 

({nX + Y},X,Y,Z) =4 (U,X,Y,Z) (9) 

where U is uniform on the unit interval independent of (X, Y, Z), if and only if fx 
is Rajchman. 

Proof. If ji is a probability measure on T 1 x R m , let us put 

m0 = J e^ kx+< ^>fi(dx,dy), 

then ii n =U n iff ji n (k, () -> fr(k, () VA; e Z, V( E R m 

a) If Fx € 72 

^{nX+Y},x,Y,z)(k, Ci, C2, Cs) = E[exp{2mA;(nX + F) + <i^T + i( 2 Y + i < Cs, Z >}] 

^ k ™f{x)V {x} {dx) 

with /(x) = E[exp{2mA;r + i( ± X + i( 2 Y + i < (3, Z >}\{X} = x\. The fact that 
£ 72 gives the result. 

b) Conversely, taking (k, (1,(2,(3) = (1,0,-2tt,0) gives F {x} (n) -> i.e. P x e 
72. o 
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The preceding definitions and properties extend to T d : a measure \i on T d is 
said to be in 72, if (i(k) — > as k — > oo in Z d . The set of measures in 72. is a 
band. If X is Revalued, limi^oo Ee l<n,x> = is equivalent to P{x} £ 72 where 
{x} = ({xi},..., 

III. Finite dimensional cases. 

In this part we carry on deeper with the basic example (section 1.2) in the finite 
dimensional case. 

1 1 1.1. We suppose Y is Revalued, measured with an equidistant graduation cor- 
responding to an orthonormal rectilinear coordinate system, and estimated to the 
nearest graduation component by component. Thus we put 

Y n = Y + -6(nY) 
n 

with%) = (i-{ yi },-..,i-{y4). 

Theorem 5. a) Iffy is Rajchman and if X is M™ -valued 

(X,n(Y n -Y)) =^ (X,(V u ...,V d )) (10) 

where the Vi 's are independent identically distributed uniformly distributed on (— |, |) 
and independent of X . 
For all up E C 1 n lip(R d ) 

(XMv(Y n )- V (Y))) =4 (X,J2VM(Y)) (11) 

i=i 

n 2 E[(v,(r„) - <p(Y)) 2 \Y = y\ -> 1 £ m ( 12 ) 

i=l 

m particular 

n*E[( V (Y n ) - V {Y)f] - E y [^ J>f (y)]. (13) 

i=i 

zs o/ c/ass C 2 ; i/ie conditional expectation n 2 K[<p(Y n ) — tp(Y)\Y = y] 
possesses a version n 2 (ip(y+ -6{ny)) — (p(y)) independent of the probability measure 
P which converges in the sense of distributions to the function ^ A ip. 
c) //Py < on M d ; G ^([0, 1]) 

(x,^(n(y B -y))) =^(x,^(v)). (14) 
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d) We consider the bias operators on the algebra C 2 of bounded functions with 
bounded derivatives up to order 2 with the sequence a n = n 2 . //Py G 7£ and if one 
of the following condition is fulfilled 

i) Mi = l,...,d the partial derivative <9jPy in the sense of distributions is a 
measure <C Py of the form p«Py with pi e L 2 (F Y ) 

ii) Py = /il G |j with G open set, h £ H l fl L°°{G), h>0 
then hypotheses (HI) to (H4) are satisfied and 

A[<p] =^A^ 

AM = Ta A V + h E ViPi case i) 

A M = M Ai P + M^ h 'i l Pi case ii) 

m = ^Evf- 



Proof. The argument for relation (10) is similar to one dimensional case. The 
relation (11) comes from the Taylor expansion tp(Y n ) — <p(Y) = 

= Ei=l(^n,i ~~ ^i) Jo fi0^ri,l, ■ ■ ■ , Yn,i-1, ^ + t(Y n j — Yj), F+l, . . . , Yd) dt 

and the convergence 



(X,J20(nYM(Y)) 



thanks to (10) and the following approximation in L 1 



E 



Jo 



)dt 



To prove the formulas (12) and (13) let us remark that 
n 2 E[^(Y n ) - V {Yf\Y = y] = 



= E 



V OinYi) f . . , Yi + t(Y n , - Yi), . . )dt 
Jo 



\Y = y 



^ / l I 

y^0{nyi) \ y' i (y 1 + -0(ny 1 ),...,y i + t-0(ny i ),...)dt 
^ In n n 



0. 



Py - a.s. 



each term (9(nyi) </^(. . .)dt) 2 converges to J 9 2 (p' 2 (y) = j^tp' 2 in L 1 and each term 

9(nyi)9(nyj) Jq • • • Jq ■ ■ ■ goes to zero in L 1 what proves the a) of the statement. 

The point b) is obtained following the same lines with a Taylor expansion up to 
second order and an integration by part thanks to the fact that (p is now supposed 
tobeC 2 . 



11 



In order to prove c) let us suppose first that Py = l[ 01 ]d.dy. Considering a 
sequence of functions ipk £ Cb tending to ip in L 1 we have the bound 

|jgj e i<u,X> e ivip(6(nY))^ _ ^ e i<u,X> e iviP k (d(nY))^ 

< \v\J\ip(e(ny))-M^ny))\dy 

= M e;= • • • IT 1 ■ ■ ■ \m*yi) ■■■)- M0(nVi) ■ ■ -)\dyi ...dy d 
HHE-E/''JI*),---)-#(4---)l^-^ 



This yields (14) in this case. Now if Py <C dy then P{y} <C dy on [0, l] d and the 
weak convergence under dy on [0, l] d implies the weak convergence under P{y} what 
yields the result. 

In d) the point i) is proved by the approach already used in Bouleau (2006) 
consisting of proving that hypothesis (H3) is fulfilled by displaying the operator A 
thanks to an integration by parts. The point ii) is an application of Girsanov theo- 
rem (theorem 4). o 

Remarks. 1) About the relations (11) (12) (13), let us note that with respect to 
the form 

i 

when it is closable, the random variable ^ appears to be a gradient : if we 
put y? # = J2i ViV'i then a we have 
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the square field operator associated to £. We will find this phenomenon again on 
the Wiener space. 

2) If d — 1, the convergence in (12) holds in L p 1 < p < oo and in d) belonging 
to 1Z is automatic under i) or ii). 

3) Let us also remark that when d = 1 assumptions i) and ii) may be replaced 
by the Hamza condition on Py which suffices to imply (HI) to (H4). 

1 1 1.2. Approximation to the nearest graduation, by excess, or by default. 

Let us come back to the basic example. When the approximation is done to 
the nearest graduation, on the algebra C% the four bias operators are zero with the 
sequence a n = n (cf. theorem 5 where a n = n 2 ). 

But we would obtain a quite different result with an approximation by default 
or by excess where the effect of the shift is dominating. 

If the random variable Y is approximated by default by = then 



n(YW-Y)=U-U and - Y)} -> -\ 
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as soon as Y is say bounded. With this approximation, if we do not erase the shift 
down proportional to — and if we take a n = nwe obtain first order bias operators 
without diffusion : A[(p] = —\ip' = — A[(p] and A = 0. The same happens of course 
with the approximation by excess. 

1 1 1.3. Extension to more general graduations. 

Let Y be an Revalued random variable approximated by Y n = Y + £ n (Y) with a 
sequence a n \ oo on the algebra V = C{e <u ' x> , u G W 1 }, the function £ n satisfying 

f a n E[\U 3 (Y)] - 



0< 



a n E[ip(Y) < u,£ n (Y) > 2 ] -»• E Y [ip.u*ju] Vy? G T>, \/u G 

-2 

■ 111 U1SW1UUUU11S SC11&C c ± 



with 7y G L°°(P y ) and |]Sf in distributions sense G L 2 (Py) 



k aj%(y) < u,UY) >] ^ V</? G I?. 
Under these hypotheses we have 

Theorem 6. a) (HI) is satisfied and 

Tr i 1 V- 9 V 



6j If for i = 1,. . . ,d, the partial derivative <9jPy in the sense of distributions is a 
bounded measure of the form p^Py with pi G L 2 (Py) then assumptions (HI) to (H4) 
are fulfilled and Wip G D 



dxidxj 



<9<p 



t/ie square field operator is 



m = E 



<9<p <9<p 



Proof. By the choice of the algebra X> this theorem is simple. The argument con- 
sists of elementary Taylor expansions in order to prove the existence of the bias 
operators. Then theorem 2 applies. o 

IV. Rajchman martingales. 

Let {!F t ) be a right continuous filtration on (f2, A, P) and M be a continuous local 
(jF t , P)-martingale nought at zero. M will be said to be Rajchman if the measure 
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d(M, M) s restricted to compact intervals belongs to 1Z almost surely. 

IV. 1. We will show that the method followed by Rootzen (1980) extends to Rajch- 
man martingales and provides the following 

Theorem 7. Let M be a continuous local martingale which is Rajchman and s.t. 
(M, M)oo = oo. 

Let f be a bounded Riemann-integrable periodic function with unit period on H. 
s.t. /q f(s)ds = 0. We put T n (t) = inf{s : f°f 2 (nu)d(M,M) u > t}. Then for 
any random variable X 

(X, [ f(ns)dM s ) =4 (X,W limM , M) ), (15) 
Jo 

the weak convergence is understood on R x C([0, 1]) and W is an independent stan- 
dard Brownian motion. 



Before proving the theorem, let us remark that it shows that the random measure 
dM s behaves in some sense like a Rajchman measure. Indeed if Py G 1Z we have 

/■i rv 

g(nx)¥ Y {dx) g(x)dx F Y (dx) 

> JO J -oo 

as soon as g is periodic with unit period, Riemann-integrable and bounded. Now 
applying the theorem to the Brownian motion gives the similar relation 

f fins) dB s =4 ( T f 2 {s)dsf' 2 f dW s . 
Jo Jo Jo 

Proof. We consider the local martingale N t = J* f(ns)dM s . 

a) In order to be sure that (N, N)^ = oo, we change N t into N t = J Q * f n (s)dM s 
with f n (s) = f(ns) for s G [0, 1), / n (s) = for s G [l,n] and f n (s) — 1 for t > n. 
We put S n (t) = inf{s : (N,N) S > t}. 

b) We want to show 

E[£F(7V 5 J] - EKF(W)] V^GL^P) VF G C 6 ([0, 1]). (16) 

It is enough to consider the case £ > 0, E£ = 1, and £ may be supposed to be 
jF T -measurable for a deterministic time T large enough. Let be P = £.P and -D(t) = 
E[£|jF t ]. The process 



M t = M t - / D-\s)d(M,D c ) s 
Jo 



is a continuous local martingale under P. Therefore Jq 5 ™^ c^^s is a Brownian 
motion under P Revuz and Yor (1994, p. 313 theorem 1.4 and p 173). Writing 

fS n (t) rS n (t) rS n (t) f l \ 

\ f n {s)dM s = f n (s)dM s + ^->d(M,D c } s 

Jo Jo Jo L'K 8 ) 
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and noting that d(M,D c ) s vanishes on ]T, oo[, in order to show (16) it suffices to 
show 

* US) 



sup 

0<t<T 



hence to show 



sup 

0<t<l 



D(s) 

1 f(ns) 



d(M, D c 



d{M,D c ) t 



a.s. when n — > oo 



a.s. when n — > oo 



/o D(s) 

and, because M is Rajchman this comes from the following lemma : 
Lemma. Let f be as in the statement of the theorem, then V/x G 1Z 



sup 

0<t<l 



f(ns)n{ds) 



as n — > oo. 



Proof. We have 



f(ns)fi(ds) 



f f(s)ds [ n(ds) = 0. 
o Jo 



Since / is bounded, the functions J* f(ns)fj,(ds) are equi-continuous and the result 
follows from Ascoli theorem. o 

c) This proves the following stable convergence 

(X, / f(ns)dM s ) =4 {X,W) 
Jo 

and by the fact that the following limit 

f f(ns)d(M, M) s ^ f f 2 (s)ds(M, M) t 
Jo Jo 

is a continuous process, this gives the announced result. o 
Remark. If J 1 f(s)ds ^ 0, then keeping the other hypotheses unchanged, we obtain 



(X, / f(ns)dM s ) 
Jo 



X,( f{s)ds)M. + { {f- ff) 1/2 W {MM) ) . 



IV. 2. Limit quadratic form for Rajchman martingales. 

We study the induced limit quadratic form when the martingale M is approxi- 
mated by the martingale M™ = M t + J* * ^f(ns)dM s . The notation is the same as in 
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the preceding section and / satisfies the same hypotheses as in theorem 7. 

Theorem 8. Let M be a Rajchman martingale s.t. Mi e L 2 and r\, ( bounded 
adapted processes. Then 



n 2 E 



(exp{i I r] s dM^}-exp{i [ r] s dM s })(exp{i [ ( s dM?} - exp{i I ( s dM s }) 
Jo Jo Jo Jo 



-E 



3xp{* I (rj s + Cs)dM s } [ r ]s ( s d(M,M) l 
Jo Jo 



1 

2 



f 2 (s)ds. 







Proof. By the finite increments formula, the first term in the statement may be 
written 

-E[exp{* [ (r] s + QdM s } [ Vs f{ns)dM s [ ( s f(ns)dM s ] + o(l) 
Jo Jo Jo 

therefore, thanks to theorem 7, the theorem is a consequence of the following lemma : 

Lemma. Suppose EMf < oo and r\ adapted and bounded, then the random variables 
Jo Vsf( ns )dM s are uniformly integrable. 

Proof. It suffices to remark that their L -norm is equal to E f* r] 2 s p{ns) d(M, M) s 
hence uniformly bounded. o 

V. Sufficient closability conditions on the Wiener space. 

The closability problem of the limit quadratic forms obtained in the preceding 
section, may be tackled with the tools available on the Wiener space. 

Let us approximate the Brownian motion (B t )te[o,i] by the process B™ = B t + 
Jo nf( ns ) dB s where / satisfies the same hypotheses as before. We consider here 
only deterministic integrands. 

Theorem 9. a) Let £ e L 2 ([0, 1]), and let X be a random variable defined on the 
Wiener space, i.e. a Wiener functional, then 

X,n(exp{t [\dB n }-exp{t f\dB})) =4 (x, ||/|U 2 (exp{i [\dB})# 

v Jo Jo J \ Jo 

(17) 

here for any regular Wiener functional Z we put Z # (u,w) = j Q D S Z dW s , where W 
is an independent Brownian motion. 

b) 

n 2 E [ (e *B» _ e ^.B )2 ] ^ _ E [ e 2i€.B] f\ 2 ds\\f\\ 2 L2 (18) 

Jo 
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on the algebra C{e^ B } the quadratic form — |E[e 2 ^' B ] Jq $, 2 ds is closable, its closure 
is the Ornstein- Uhlenbeck form. 

Proof, a) The first assertion comes easily from the similar result concerning Ra- 
jchman martingales using the fact that e ta ^ n^ ns ) dBs da — > 1 in L' p p G [1, oo[. 

b) The obtained quadratic form is immediately recognized as the Ornstein- 
Uhlenbeck form which is closed. It follows that hypothesis (H3) is fulfilled and 
the symmetric bias operator is 

A[e^n = (-*- jtidB-\j eds^j 

o 



If instead of the Wiener measure m, we consider the measure mi = h.m for 
an h > 0, h G D ou H L°° where 3 ou (= D 2,1 ) denotes the domain of the Ornstein- 
Uhlenbeck form, we know by Girsanov theorem (theorem 4) that the form — ^Ki[e 2l ^ B J* £, 2 ds] 
is closable, admits the same square field operator on D ou , and that its generator A 1 
satisfies 

A^ip] = A[(p] + 2}fcm[p, h] for y? G VA ou 

Since the point a) of the theorem is still valid under m\ because of the properties of 
stable convergence, the preceding theorem is still valid under mi, the Dirichlet form 
being now 

S 1 [^] = ^E 1 [T ou [ V ]] for y?GD ou . 

Remark. Let us come back to the general case of Rajchman martingales. If we 
suppose the Rajchman local martingale M is in addition Gaussian, which is equiva- 
lent to suppose (M, M) deterministic, then on the algebra C{e l ^ dM ; £ deterministic 
bounded } the limit quadratic form 

_ E[e i/(^K)«w f 1 ( s r] s d(M, M) s ] H/IH2 
Jo 

is closable, hence (H3) is satisfied. 

Indeed, it suffices to exhibit the corresponding symmetric bias operator. But 
by the use of the calculus for Gaussian variables, it is easily seen that the operator 
defined by 

Aft**™] = j^ dM J^dM-^J ed(M,M))j j f 2 ds 
satisfies the required condition. o 
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In the case f(x) = 6{x) — | — {x}, the approximation used in the parts IV 
and V consists in approximating B t by B t + J"* ^9(ns)dB s . It is the most natural 
approximation suggested by the Rajchman property and the arbitrary functions 
principle. It yields also other approximation operators on the Wiener space, cf. 
Bouleau (CRAS 2006). 

But it is different from the approximations usually encountered in the discretiza- 
tion of stochastic differential equations. 

In order to draw a link between the preceding study and works concerning the 
discretization of stochastic differential equations by the Euler scheme, especially 
those of Kurtz and Protter (1991), and Jacod and Protter (1998), we may remark 
that the preceding results which yield 

(n J\s - [ ^)dB s , n J\B a - B^)ds, B^j =4 {-±=W.+\b, --±=W.+\b, B) 

(19) 

are generally hidden by a dominating phenomenon 

(v^J o (B s -B^)dB s ,B^ =U(±W.,,B.) 

due to the fact that when a variable of the second chaos (or in further chaos) 
converges stably to a Gaussian limit, this one appears to be independent of the first 
chaos and therefore of B itself. 

The stable convergence (19) acts even on the first chaos. It concerns, for example, 
stochastic differential equations of the form 

f X\ = 4 + JJ r{Xl)dB s + JJ r (Xj, X 2 )ds 

I x^xl + fj^xlxDds [ Uj 

where X 1 is with values in IR fcl , X 2 in R fe2 , B in M. d and / u are matrices with 
suitable dimensions. Such equations are encountered to describe the movement of 
mechanical systems under the action of forces with a random noise, when the noisy 
forces depend only on the position of the system and the time. Typically 

r x t = x + j^v s ds 

\V t = V + Jj a(X s , V 8 , s)ds + b(X s , s)dB s 

which is a perturbation of the equation = a(x,^,t). In such equations the 
stochastic integral may be understood as Ito as well as Stratonovitch. For the equa- 
tion (20) the method of Kurtz and Protter (1991) without major changes yields the 
following result that we state in the case k\ — k 2 — d — 1 for simplicity. 

Theorem 10. If functions / u are C\, and if X n is the solution of (22) by the Euler 
scheme, 

(n(X n -X),X,B)^(U,X,B) 
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where the process U is solution of the stochastic differential equation 

u v = E f d -&^ u " dY ' - E f §r« E f km WdZT j 

k,j J ° k k,j J ° k m 

where Y s = (B s , s)' and 

dz s' = -^k dw s + ^ 
dzf = f 

and as ewer IF zs an independent Brownian motion. 
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